set scheme s1color 

use TeO, clear

// Standardize outcomes
foreach x in i_ecole harass_school {
	qui sum `x'
	local `x'm=r(mean)
	sum `x' 
	local sd=r(sd)
	gen `x'st=(`x'-``x'm')/`sd'  
}


//Regressions
local sample if bpl=="1101" &female==1  
local inter interm 
local controls i.birthyear i.religion tmuslim 
estimates clear

reg harass_schoolst `controls' `inter' `sample', cl(religion)
eststo m1

reg i_ecolest `controls' `inter' `sample', cl(religion)
eststo m2

// Graph
coefplot (m1, aseq("Experienced racism in school") ciopts(lcolor(black) recast(rcap))) ///
(m2, aseq("Trust in French school") ciopts(lcolor(black) recast(rcap))), ci(90) xline(0) keep(interm) ///
plotregion(style(none)) ysca(titlegap(2)) msize(small) mcolor(black) ///
xtitle("Estimated difference in differences", size(small)) ytitle("", size(small)) ///
xlabel(, nogrid labsize(small)) ylabel(, nogrid labsize(small)) xsca(titlegap(2)) ///
swapnames legend(off) recast(bar) barwidth(0.25) color(blue) fcolor(*.5) citop
